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Abstract 

The diamagnetism of confined Dirac fermions submitted to a uniform magnetic field in dis- 
ordered graphene is investigated. The solutions of the energy spectrum are used to discuss the 
orbital magnetism from a statistical mechanical point of view. More precisely, by the technique 
of Green functions the self-energy for short and long-ranged disorders is obtained. This allows us 
to determine the susceptibility for short and long-ranged disorders together with confinement. We 
compare our results with already published work and point out the relevance of these findings to a 
systematic formulation of the diamagnetism in a confining potential. 
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1 Introduction 



The magnetism of graphene was first studied as a simple model for three-dimensional graphite [T] where 
the susceptibility of the disorder- free graphene was calculated within the effective mass approximation. 
It was found that the system exhibits a large diamagnetism at the Fermi energy ep = 0, expressed as 
a 5- function of ep at the absolute zero temperature. The graphene magnetism was considered again 
in studies on the graphite intercalation compounds, where the tight-binding model was applied for a 
wide range of Fermi energies [21 |3l [U [5] . The effects of disorder on graphene under magnetic fields have 
been examined in early theoretical studies before the experimental discovery of graphene, where the 
electronic structure [B] , the transport properties [3 El [B] , and the de Haas- van Alphen effect [S] were 
investigated. Recently the Shubnikov-de Haas oscillation was studied in disordered graphene [TOl [11] 
and the spectral and transport properties were examined in presence of lattice defects under the 
magnetic fields [12] . 

Very recently, an interesting development on the diamagnetism of disordered graphene was reported 
by Koshino and Ando [13j. They studied the graphene monolayer orbital magnetism within the 
effective mass approximation. In models of short and long-ranged disorders, the magnetization was 
calculated with the self-consistent Born approximation. In the zero-field limit, the susceptibility 
becomes highly diamagnetic around zero energy, while it has a long tail proportional to the inverse 
of the Fermi energy. It was demonstrated how the magnetic oscillation vanishes and converges to the 
susceptibility, on going from a strong- field regime to zero field. Additionally, the behavior at zero 
energy was shown to be highly singular. 

On the other hand, an exact solution of a related problem, that has been studied at various levels 
by researchers dealing with different physical issues (see for instance [14:\ [T5] ) , was given by one of us 
(A J) and his collaborators [16], considering a relativistic particle subjected to an external magnetic 
field as well as to a confining potential. By a similarity transformation the system can be diagonalized 
in a simple way. Solving the eigenvalue equation, and accounting for the complete space of the 
eigenfunctions one can include various cases related to different physical settings. More precisely, 
from the nature of the problem it was possible to obtain separate angular and radial solutions. The 
radial equation lead to an exact relation between the two-spinor components. In fact, depending on 
the range of values of three physical quantities, the full solution space split into eight disconnected 
subspaces as summarized in table 1. 

Motivated by different investigations on Dirac fermions in (2 -|- l)-dimensions, in particular by 
references \13\ 116], we treat diamagnetism of a confined system in a statistical mechanical way. More 
precisely, we study the orbital magnetism of Dirac fermions in uniform magnetic field, disordered 
graphene and confining potential. In fact, we combine studies reported in [13\ [T6] to generalize the 
results of Koshino and Ando [13] about diamagnetism in disordered graphene to the confinement case. 
This can be done by using the energy spectrum solutions to study the self-energy for both regimes: 
short and long-ranged disorders by the technique of Green functions. The self-energy allows us to 
obtain the shape of the density of states that is needed to determine the related thermodynamical 
quantities and discuss different issues. 

Subsequently, we calculate the susceptibility for two regimes and underline what makes the dif- 
ference with respect to the standard case [13], namely the analysis without confinement (k = 0). 



1 



For short-ranged disorder, we obtain interesting results in terms of a parameter of confinement and 
disorder strength, called C . In particular, we show that there is a quantum correction to the result 
obtained by Koshino and Ando [13], which disappears by switching off k. Furthermore, we notice that 
the clean limit result without confinement can be obtained by considering C going to zero. 

As far as the long-range disorder is concerned, the susceptibility shows an additional second term 
of the order of 0{C) compared to the short-range case. However, this gives a minor effect since C is 
assumed to be small. When the terms of the order 0{C)'^ are neglected, the susceptibility becomes just 

1 — 3C times as large as in the short-ranged disorder. Accordingly the integration of the susceptibility 
over the energy e depends weakly on C, while in the limit C — )• we again get the susceptibility 
as (5-function. Finally, we notice that the case k = allows us to recover the results and related 
conclusions proposed in jl3j . 

The paper is organized as follows. In section 2, we review the energy spectrum solutions of the 
confinement problem needed to deal with different issues. Section 3 is devoted to introduce disordered 
graphene where we give the corresponding four component spinors as well as the disorder and confining 
potentials. In section 4, we use the self-consistent Born approximation to determine the self-energy 
and the density of states. These will allow us to treat the orbital magnetism by distinguishing between 
the short- and long-ranged disorders in section 5. Finally, we conclude by discussing the main results 
and possible extension of our work. 

2 Confinement problem 

We start by formulating the problem in terms of our apparoach. In doing so, we introduce a similarity 
transformation of the Dirac equation in polar coordinates. This will be convenient to handle the exact 
relationship between spinor components and thus derive the full spectrum accounting for the complete 
Hilbert space. 

2.1 Hamiltonian 

The problem of a charged particle moving in a constant magnetic field B = B z is a. 2D problem in 
the plane normal to the field [the Cartesian (x,y)-plane or cylindrical (r, ^)-plane]. In the relativistic 
units, h = c = 1, the Dirac equation in (2 + l)-dimensions for a spinor of charge e and massless in the 
electromagnetic potential = {Aq,A) reads as follows 

[il^{d^ + ieA^)]ij = 0, ^ = 0,1,2 (1) 

where the summation convention over repeated indices is used. 7^ = (7'', 7) are three unimodular 
square matrices satisfying the anti-commutation relation: 

i^^^^^Y}=rY + l''l^ = 2g^'^ (2) 

where Q is the metric of Minkowski space-time, which is equal to diag(+ — — ). A minimal irreducible 
matrix representation that satisfies this relation is given by 7^ = (T3, 7 = i (j where are the 

2x2 hermitian Pauli spin matrices: 
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Equation ([T]) can be rewritten as 



d 

i 11^=1 -ia-V + ea-A + eAo\'il} (4) 



dt 

where a is the hermitian matrix a = i cr^a. We will see below that the symmetry of the problem 
is preserved even if we introduce an additional coupling to the 2D Dirac-oscillator potential. This 
coupling is introduced by the substitution V — )• V + Kra^ where k is a constant parameter. 

For time independent potentials, the two-component spinor wavefunction r, 9) is written as 

i;{t,r,e) = e-''*ij{r,e) (5) 

and dH) becomes the energy eigenvalue equation (T-L — e)^ = where e is the relativistic energy. The 
Dirac Hamiltonian T-L is the 2x2 matrix operator 

H = ia^a ■ flir + i<J3a ■ OTie (6) 
where (f, 9) are the unit vectors in cylindrical coordinates and 

Tir = —idr + eAr — zKras, (7) 

Ue = — de + eAg. 
r 

For regular solutions of (HD, square integrability (with respect to the measure d?f = r dr dO) and the 
boundary conditions require that il:[r,9) satisfies 

V^i^{r,0)l=^ = Q, y/^^{r,9)l_^^ = Q, ^(0 + 27r) = ^. (8) 

To simplify the construction of the solution, we look for a local 2x2 similarity transformation 
A(r, 9) that maps the cylindrical projection of the Pauli matrices (<? • r, a ■ 9) into their canonical 
Cartesian representation (o"i,(72), respectively [T7]. That means 

Aa-fA"^=o-i, Aa-eA~^=a2. (9) 

A 2 X 2 matrix that satisfies this requirement is 

A(r,^) = A(r,^)e5'^3e (^q) 

where A(r, 9) is a real function and the exponential is a 2 x 2 unitary matrix. The Dirac Hamiltonian 
([6]) gets mapped into 

H = A^A-i = -a2Hr + aiHe (11) 



where 



Hr = —i — j-^ + ieAr — iKra-s, (12) 
He = -Ude-^- ^CJ3 j + eAe 



with Afc = 9fcA. Therefore, the 2x2 Dirac Hamiltonian becomes 



H = 



<9r - X + 27 + ^eA- ~Kr~^{de-^)+ eAg 

-dr + IT - ^ - i^^r - i^r- ^{de- ^) + eAe 



(13) 



3 



Thus, hermiticity of (|13p requires that 

A.=0, ^-1=0 (14) 

and fixes the exact form of the modulus of the similarity transformation to be X{r,9) = y^. It is 
interesting to note that turns out to be the integration measure in 2D cylindrical coordinates. We 
could have eliminated the A factor in the definition of A in ()10p by proposing that the new spinor 
wavefunction x be replaced by -^x(^)^)- that case, the transformation matrix A becomes simply 
e2'^^^ , which is unitary. However, making the presentation as above gave us a good opportunity to 
show (in a different approach) why it is customary to take the radial component of the wavefunction 
in 2D cylindrical coordinates to be proportional to Finally, we obtain the (2 + l)-dimensional 
Dirac equation (H — e) x = ^ for a charged spinor in static electromagnetic potential as 



-£ dr + ieAr — Kr — ^de + eAe 



-dr - ieAr - Kr - ^de + eAe 




(15) 



where x± ^-re the components of the transformed wavefunction |x) = A\iIj). This equation will be 
solved by choosing an appropriate gauge to end up with the full Hilbert space. 

2.2 Energy spectrum 

Now, we specialize to the case where a constant magnetic field of strength B is applied at right angles 
to the (r, 0)-plane, which is B = B z. Therefore, the electromagnetic potential has the time and space 
components 

Ao = 0, A{r,e) = ^Br9. (16) 
Consequently, (jlSp becomes completely separable and we can write the spinor wavefunction as 

X±{r,e) = Mr)T{9). (17) 

Thus, the angular component satisfies —i^ = where ^ is a real separation constant giving the 
function 

.W = ^e.^^ (IS) 



On the other hand, the boundary condition tl'{9 + 27r) = ^(0) requires that e*2<e-''^3'' = +1 which, 
in turn, demands that e*^'^^ = — 1 giving the following quantum number: 

^ = ^v, v = ±l,±3,±5--- . (19) 
Consequently, the Dirac equation for the two-component radial spinor is reduced to 



where the physical constant u is given by w = Ir?' — k and Ib is the magnetic length Ib = 
Thus, the presence of the 2D Dirac-oscillator coupling did, in fact, maintain the symmetry of the 
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problem as stated below (jH). Moreover, its introduction is equivalent to changing the magnetic field 
as eB — > eB — 2k. As a result of the wave equation (j20p . the two spinor components satisfy the 
"kinetic balance" relation 

Mr) (21) 



h{r) = - 



T-r "I 

ar r 



£ 

where e 7^ 0. Therefore, the solution of the problem with the top/bottom sign corresponds to the 
positive/negative energy solution. Using the exact relation (|2ip to eliminate one component in terms 
of the other in ()20p results in the following Schrodinger-like differential equation for each spinor 
component: 

{-^ + + + [-e- + ui2^± D] } 0i(r) = 0. (22) 

We stress that this equation gives only one radial spinor component. One must choose either the top 
or bottom sign to obtain the component that corresponds to the positive or negative energy solutions, 
respectively. The second component is obtained by substituting this into the exact relation ()2ip . 
Nonetheless, we only need to find one solution (the positive- or negative-energy solution), because the 
other is obtained by a simple map. For example, the following map takes the positive energy solution 
into the negative energy solution: 

e — > —e, V — > —V, bj — > —uj, (j)± — > (/)=p (23) 

which, in fact, is the CVT transformation. Here the charge conjugation C means that e — > — e and 
K — 7- —K. It is easy to check that the above map (|23p originates from the fact that the Dirac equation 
(|20p is invariant under such transformation. Hence, we just need to solve for positive energies and 
use the above transformation to obtain the negative energy solutions. The total spinor wavefunction 
reads 

Mr, e) = ^ e*«^e-i'^»V(^) (24) 
where (j){r) has two components, such as 

Equation (|22p looks like the non-relativistic oscillator problem with a certain parameter map of the 
frequency, angular momentum, and energy. For regular solutions of (|22p . the bound states will be of 
the form 

</.± ~ z'^e-^/2L^(z) (26) 

where Ll^{z) is the associated Laguerre polynomials of order n = 0, 1,2, ••• and z = p^r"^. The 
constants {/i, v, p} are real and related to the physical parameters B, k, and ^. Square integrability 
and the boundary conditions require that 2p > ^ and u > —1. 

Substituting the ansatz (|26p into ()22p and using the differential equation for the Laguerre poly- 
nomials |18j . we obtain four equations. Three of them determine the parameters {/i,z^, />} and one 
determines the energy spectrum. The first three are 

1 

2n 

z 

( f: - I ^ n 

(27) 
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For regular solutions of ()22p . the it sign in the expression for v corresponds to it^ > 0. The fourth 
equation gives the following (positive and negative) energy spectra: 



2n + 1 ± ^-^ + i{s + s>) 



(28) 



where s = sgn (w) = and s' = sgn {^). The sign of uj depends on whether k is larger or smaller 
than the magnetic length Ib- To compare our work with frequently used notation in the literature, 
we can replace the quantum number ^ by A; + ^, where fc = 0, ±1, ±2, • • • and ^ — > — ^ implies that 
k — 7- —k — 1. In that case, one may write the positive eigenvalues as 



e 



n,k 



^/2\oj\ [2n + l + s + k{s + s')] (29) 



and the negative ones as 



e 



n,k 



-yj2\oj\[2n + I + s' + k{s + s')] (30) 



where s' = +1 for /c = 0. It is interesting to note that for < the spectrum is infinitely degenerate 
because it is independent of ^. However, for .^o; > the degeneracy is finite and equal to n + k + 1. 
Substituting the wavefunction parameters given by (|27p into the ansatz (j26p gives for e > 0: 



as well as for e < 0: 



, ^1^2 { AttLiix^) , k>0 

I K^k Ln'~' i^') , k<0 

where x = ry^M" are normalization constants that depend on the physical quantities Ib 

and K. The lower component (j)-{r) is obtained by substituting (|3ip and (|32p into the exact relation 
(|2ip . Doing so while exploiting the differential and recursion properties of the Laguerre polynomials 
we obtain the following for e > 0: 



(j)_(r) = -^^q^a;l'^+2l e 2^ 



<1 X [(« - l)^n {x') + 2L^i (x^)] , k>0 

k < 0. ^^^^ 



^n,k 



{s - l)(n - fc)L-^-i (x2) - (. + l)(n + l)L-^ri (x' 



On the other hand, repeating the same calculation for the upper component of the negative energy 
solution gives the function 



(j)^{r) = -vi^^^l xl^+sl e~^^^ 



^n,k 



K% [(1 + s)(n + /c + \)Li (x2) + (1 - s)(n + \)Li^^ (x2)] , /c > 
A;^^x[(1 + .)L-^-1(x2)-2L-^^(x2)], fc<0 



which can have also been obtained by applying the CVT map (j23p to (j33p . Thus, the structure of 
the whole Hilbert space solution consists of eight disconnected spaces that can be displayed in tabular 
form as shown in Table 1: 
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Frequency 






Energy 




e < 


£ > 


e < 


Azimuth 


k>0 


k<0 


k>0 


k<0 


A; > 


A; < 


> 


A; < 



Table 1: Complete space solution. 

Using the standard definition, we calculate all normalization constants in the above wavefunctions. 
These are summarized in the Table 2, where as stated above s = sgn(a;) = |a;|/a;: 



Energy 


Azimuth 


Nor malizat ion 


e > 


A; > 


Km - y 


' 2n! 
7r(n+A;)! 


^{2(n + fc + l) + n(l-.)}j ' 


e > 


A: < 


<fc = y 


' n! 
7r(n— fc)! 


-^{2(n + l) + (fe + l)(s-l)}l ' 


e < 


A; > 


Kt ^ 


/ n\ 
7r(n+fc+l)! 


-^{2(n + l)+A;(s + l)} 




e < 


A; < 


^n,fc = y 


/ 2n! 
7r(n— fc— 1)! 


7^{2(n-fc)+n(. + l)} 


-1 



Table 2: Normalization in terms of different physical quantities. 

3 Disordered sublattices 

Our main task is to analyze the diamagnetism of confined Dirac fermions in disordered graphene. 
In the previous section, we settled the required tools for the confinement, however we still need to 
introduce the disorder potential and related matter. For this, we write the eigenspinors as well as the 
impurity and confining potentials in four components. 

Having obtained the energy spectrum solutions corresponding to one sublattice, it is worthwhile 
to deal with a system of graphene. Such a system is composed of a honeycomb network of carbon 
atoms where a unit cell contains one atom each from a pair of sublattices, denoted by I and II. It can 
be described by a 4 x 4 matrix Hamiltonian, such as 



H 








a 








\ 




at 
























at 




v 








a 





/ 



(35) 
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where in polar coordinates the two operators a and take the form 



5r - y + — + leAr - Kr - - Ide - — ] + eAg 



-9r- + ^ - 77 — ieAr - Kr - - { do 
X 2r r 



X, 



+ eAe. 



(36) 
(37) 



Clearly, they act on one component of the eigenspinors ,F{f ,Ff ,F^' where Ff,Fj^ represent 
the envelope functions at I and II sites for point K, respectively, and Ff^ , F^ for K', with K and 
K' indexing the two valleys. 

It is convenient for our task to label the eigenspinors by a = (j, n, k) with the valley index 
j = K,K' , the Landau level index n = 0, ±1, • • • and the wavevector k. The eigenspinors for A; > 
and e > can easily be deduced from the above analysis. For K we have 



( Li {^') \ 






(38) 



as well as 



2-^ 



x[{s-l)Li (x2)+2L^i (x2)] 



V 



(39) 



for the point K' . Without loss of generality in the forthcoming analysis we only focus on the case 
where k > and e > 0. Other configurations can be recovered from the first one by making use of 
some mapping. 

To consider disordered graphene and make comparisons with already published work [13] , we intro- 
duce the disorder potential characterized by two simple models: short- and long-ranged scatterers [6]. 
The first is an on-site potential localized at a particular I or II site with a random amplitude. A 
scatterer on site I at Ri is represented as 

f 1 z^z{ \ 



zjz[* 1 





U{f) 



V 



Ui5{r- Ri) 



(40) 



and that on site II at Ru as 



/ \ 



U{r) 







1 



zl,z[. 





y ziiz[{ 1 y 



Ui5{r- Rii) 



(41) 
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where we introduced zx = e*^'^^, 



gi^'-Rx with X = I and II, and Ui = {V3a'^/2)Ui with the 



on-site energy Ui. We assume that the scatterers are equally distributed on I and II sites with density 
nj = np = nj/2 and the mean square amplitude {{uj)'^) = {{uY)"^) = uf. 

Dominant scatterers in graphene are expected to have a potential range larger than the lattice 
constant for which inter- valley scattering is much smaller than intra- valley scattering. Further, realistic 
scatterers are likely to have the range comparable to the Fermi wavelength [19\ [20l [21] . In the following, 
however, we shall assume scatterers with potential range smaller than the Fermi wavelength. The 
reason is that the results are expected [13] to remain qualitatively the same and further that actual 
calculations are practically possible. 

In this long-range model, a scatterer at R is expressed by 



C/(f) 



/ 1 

1 



\0 






1 








1/ 



Ui6{r-R). 



(42) 



We assume the scatterer density rii and the mean square amplitude uf . It was shown that the transport 
properties in the short-ranged disorder and the long-ranged one are qualitatively similar [H [T] [8] . 
To complete our model we consider the confining potential as 4 x 4 matrix as well. This is 

\ 



V{r) 



( Kr 



V 





-Kr 








Kr 










-Kr 



(43) 



which of course can easily be obtained from the former study. The above established mathematical 
tools will serve us to deal with our task. More precisely, we will see how it can be used to describe the 
diamagnetism of relativistic particles under the influence of three constraints. These are: magnetic 
field, confined and disordered potentials. 



4 Self-consistent Born approximation 

As we claimed before, one of our objectives is to generalize the results obtained in [13] to the confine- 
ment case. Actually, this can be achieved, for instance, by adopting the same method as [13], based 
on application of the self-consistent Born approximation to evaluate the self energy and therefore the 
density of states. This allows us to determine some thermodynamical quantities, in particular the 
susceptibility. 

4.1 Self-energy 

To deal with different issues, we introduce the Green function that is related to the self-energy S via 
the Dyson equation. It is 

(Gaa'ie)) = S^^,Gi{e) + (e) ^ E^^„{e){G^„^,{e)) (44) 

a" 
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where the first term contains the matrix elements of the unperturbed Green function corresponding 
to the Hamiltonian considered before, such as 

GL'(e) = = S^a'Glie). (45) 

To proceed further, we make use of the self-consistent Born approximation for our system. The 
self-energy of the disorder- averaged Green function {Ga,a') can be written as 

Saa/(e) = '^{UaaiUa2a'){Gaia2{£)) (46) 
Q102 

where the symbol (• • • ) represents the average over the impurity configurations. 

The above equation can be solved by considering the range of the disorder. Specifically, in the 
short-ranged model one can show that the self-energy and the averaged Green function are diagonal 
with respect to a [13]. Furthermore, the self-energy becomes independent of a and leads to 



{Gaa'ie)) =6^a,GUe) (47) 

where Ga{s) is given by 

G^{e) = G(e, e^) = ^ (48) 

Prom ()46p and (j48p . it is clear that one has to find the appropriate solution for E. To do so, one needs 
to introduce respective approximations. 

Having described the needed tools, let us see how they can be applied to analyze the basic features 
of the present system. According the former analysis we have two contributions to the total self-energy 
S(e) = Stot(e), namely 

Stot(e) = Sdis(e) + Sconf(e) (49) 

where Sdis and Sconf correspond to the disordered and confining potentials, respectively. In the 
forthcoming analysis, we separately determine each part. Note that, what makes a difference with 
respect to the study reported in [13] is the second contribution and therefore one can see its impact 
on such study. 

As a first step, we have to evaluate the matrix elements of different potentials. For the impurity 
potential, using the eigenspinors it is straightforward to show 



{a\Ui\a) = ^/\oj\ui5oia' (50) 

where oj is the parameter introduced in (|20p and the integration is performed over the coordinates r 
and 9. This can be used together with ()46p for the short-ranged potential, to straightforwardly obtain 

2 oo oo / \ 

""'''^'^ = ^1^1 2: E e-...-s,ot(.) ^''^ 

n=-oo fc=0 ' ^ ' 

where the cutoff function g{e) is given 

9ie) = \ < (52) 

0, otherwise. 
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As far as the confining potential is concerned, the corresponding matrix elements can be evaluated 
to end up with 



(a| lb Kr\a') 



-2s(n + A; + 1) 



(53) 



This leads to the self-energy for the confinement: 



2 "-"^ "-"^ 



n=-oo fc=0 



s{n + k + l) 



2n+(s + l)(A: + l) 



(54) 



It is clear that Sconf is strongly K-dependent, which is an expected result because of the confining 
potential expression (P3|) . 

So far, we have obtained the different contributions to the total self-energy. This can be written 

as 



oo oo 



n=-oo k=0 



1 



s{n + k+ 1) 
2n + (s + l)(A: + l) 



9{£n,k) 



e — £n,k — Stot(e) 



(55) 



This expression can be simplified by choosing s = 1 and requiring the condition n + k+1 ^ ^. In this 
case, (1^51) reduces to 



Stot(e) 



C\uj\ 



CO CO 

E E 



9{^n,k) 



(56) 



where C is a dimensionless parameter defined as 



2 ^ 

C = niUj + — 



(57) 



which depends on the set of parameters. Thus, one can fix them to derive specific results and offer 
different interpretations. 

At this stage, one can inspect the above results to underline their basic properties. One way to 
do so is to look at the case k = 0, i.e. without confining the system. This simply reduces C to the 
so-called disorder strength 

CU=o = riiuj (58) 



which leads to the self-energy 



2 °° 



ff(eri 



e-En- S(e) 



where the corresponding eigenvalues are 



£n = Ib sgn(n)V|n|. 



(59) 



(60) 



This was obtained by studying the Dirac fermions in magnetic field and disordered graphene, more 
details can be found in [13]. On the other hand, comparing (|57p and (jSSp . C can be interpreted as a 
parameter of confinement and disorder strength. 
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4.2 Density of states 

As we claimed before, the density of states is strongly needed and will play a crucial role in the 
forthcoming analysis. Specifically, it is related to different thermodynamical quantities and therefore 
allows us to determine them in an appropriate way. This statement will be clarified starting from next 
section. 

For later convenience, we consider the density of states used in jl2j by dealing with some features 
of graphene. This is 

p(£) = --y\mGJe + iO). (61) 
vr 

a 

This form can be handled by fixing different conditions. For this, we distinguish between short- and 
long-ranged disorders. Returning to our results, we have 

P(^) = --7^ lmStot(e + iO) (62) 
vr C\uj\ 

for short-ranged disorder. 

Let us recall that the transport properties in the short- and long-ranged disorders are qualitatively 
similar |13j . In the last case, the self-energy and Green function have off-diagonal matrix elements 
between {j,n,k) and (j,—n,k). Thus, we obtain 

^a,a'ii:) = SjJ'Sk,k' '^n,n'S'^(e) + (5„ _„'i;°(e) . (63) 

Splitting the self-energy into two parts 

S±t = S,t±S^'ot (64) 

we derive the positive contribution 



as well as the negative one 



where C has the same form as for the short-range case, i.e. (|57|) . These parts can be used to derive 
the density of states for the long-range case. More precisely, we obtain 

Pie) = -1 J-lm[S+t(£ + iO) + S-„t(e + ^0)] (67) 
vr o |a;| 

which reduces to that obtained in [13j by switching off the confining parameter k. With this, we finish 
the derivations of the tools needed to tackle different issues. In fact, we will see how the above results 
can be applied to deal with the diamagnetism of the present system and emphasize what makes the 
difference with respect to the case without confinement. 
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5 Thermodynamic properties 



Now we show the relevance of the above tools. We focus on the study of the diamagnetism and 
proceed in the standard way evaluating different physical quantities. More precisely, we determine 
the magnetization, number of fermions, and susceptibility to describe the physical properties of the 
system. 

5.1 Thermodynamic quantities 

We recall useful definitions of different thermodynamic quantities. The magnetization is given by 

where i}{T, fi, B) is the thermodynamic potential and /i is the chemical potential. 

To determine the number of fermions we can use one of two methods. The first one is based on 
the definition 

--(f). 



to obtain the Maxwell relation 

'dM\ fdM 



- 



On the other hand, in terms of the density of states p, we have 



(70) 



M= / p{e,B)f{e)de, (71) 



where the fermionic distribution is f{e) = 1/ (l + e^^~^^^^^'^^ . 

From the above formulas, one can establish an interesting relation. After a straightforward calcu- 
lation, we get 



£ 



M = I defie) I de'^P^^ (72) 



oo 



in terms of the density of states. Furthermore, one can also obtain the magnetic susceptibility 

dM 



X 



dB 



(73) 

B=0 



The above quantities will be simplified much more by considering the self-consistent Born approxima- 
tion and fixing the type of disorder. 

5.2 Short-ranged disorder 

To calculate different quantities we specify the nature of disorder. According to (f62|) and ([72]) . the 
susceptibility in the self-consistent Born approximation for the short-ranged disorder can be written 
as 

12 r . fi ^ r . : \ g^£tot(£^g) 

X = --— I defie) I de Im . (74) 

_B=0 
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It is convenient to introduce 



X = e - Et 



which ahows us to rewrite 



/^l I OO OO f 



n=-oo k=0 



(75) 



(76) 



According to the above equations second derivatives are needed. The first-order derivative of Stot 
with respect to B gives 

1 -1 



9Stot(e,S) 



dB 



1 



d^tot{X, B) 
dX 



dJ:t ot{X,B) 
dB 



(77) 



which leads to the second as 



tot 



dB'- 



9S 



tot 



dX 



tot 2 ^^tot / C^Etot \ ^ ^tot / t-'^tot 



dB^ dXdB \ dB 



I dB 



(78) 



One can expand (|76p into a series in terms of a function h: 



E,or{X,B) = ^At 



^ OO OO 

9^0)+ E E^(2[n + fc + l]At) 



?i=— 00,7^0 fc=l 



(79) 



where Ai = Iwl = — k\ and the function has the form 



h{t) 



2Xg{Vi) 
X^-t' 



(80) 



Note that, taking k = we recover the results obtained in [13]. This tells us that those results have 
been generalized to the present case and we will see how they can be interpreted. 

To go further, we introduce some relevant assumptions. If the condition lm(X) ^ -^/M^ is valid, 
one can simplify ([79]) to 



CM 



OO OO 



n=— 00,^0 fc=l 

(At)2 



C 1^ h{t)dt -C 



where 6{At)^ is neglected. This leads to 



h'iO) + h'ioo) 



^tot{X,B) - Stot(^,0) = -^h'{0){Atf 



(81) 



(82) 



which can be used to calculate the above derivatives. Otherwise, after neglecting all terms containing 
a power of X larger than 3, we obtain 



5Stot 




1 - 


SEtot" 


-1 


dB 


B=0 


dX 




a^Stot 




1 - 


SEtot" 


-1 


dB^ 


B=0 


dX 





C 1 



C 2j_ 



(83) 
(84) 
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With the help of ([7711781) and (f83ll8il) . the susceptibihty becomes 



X 



3ttk 



/oo 1'^ ( 

defie) Im / de' il 
-oo J —oo \ 



5S 



tot 



dX' 



1 



B=0 



By introducing 



it is not hard to find 



de' 



dX' 



X 



def{e) Im- 



dX' 



1 



B=0 



[e - Stot(e)]' 

As far as the magnetization is concerned for the short-ranged disorder, one can use ()62p . ([8^ 
and ([5U|) . These give 



(85) 

(86) 
(87) 



M 



e 



def{e) Im- 



[e - Stot(e)]' 

On the other hand, the number of fermions can be also formulated as 



1 2 r°° 



(89) 



Now we will see how the above results will be simplified by making use of different considerations. 
In fact, to underline what makes a difference with respect to other studies, we consider the zero- field 
limit and determine the corresponding susceptibility. More precisely, in such limit (|56|) can be written 
as 



= C 



1 



(e-Stot)2-t2 12 (e-S 



+ 



tot J 



which can be approximated by assuming that £ <^ £c- In this situation, we end up with 



Stot(e) = -C(e-Stot)log 



{e - Stot)2 



1 



^2 

12 ^(e - Stot)3 



H . 



(90) 



(91) 



To further simplify the above form, one may consider the condition k ^ 1 and choose an appropriate 
branch of the logarithm. Thus, one gets 

-1 



Stot(e) =£-£ 



2Ch 



ie \ 

2rc) 



(92) 



where /l(-z) is the Lambert C— function (called also Omega function which is the inverse function of 
f{z) = ze^ [22]) and T is given by 

(93) 



r = .cexp (-^ 

In the region where the energy fulfills the constraint |e| ^ F, Stot becomes 



Stot(e + «0) « -2eClog 



i-K \e\ C. 



(94) 



Following from this expression, we can derive specific results and offer different discussions. Let us 
split dMl) as 



Stot(e + «0) 



2e log 








2elog 








£ 




4 




£ 





(95) 
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The first term is similar to that obtained in [13] and becomes exactly the same if the constraint k = 
is taken into account. Accordingly, we can interpret the second as a quantum correction. 

To treat the susceptibility for short-ranged disorder, we distinguish two cases. In the first case the 
energy is constrained by the condition e ^ Ec- At zero temperature, we evaluate the integral to end 
up with 



=2 



where the function F{x) is given by 

1 



fU-ix) + ^fli-ix) 



(97) 



F(x) = Im 

X 

Without confinement, this result reduces to 

which is similar to that obtained in J3j. Therefore, ()96p is general in the sense that one can change 
two parameters, i.e. disorder Ui and confinement k, to offer different interpretations. In particular, 
we mention that even when the disorder becomes smaller, the peak of (j96p does not become narrower. 
However, it happens when C becomes smaller. On the other hand, by noticing that F{x) has its 
maximum at x = with F{0) = 1, we obtain 

p2 or 



It is worthwhile to see what happens in the case where \e\ ^ F. Using ()94p the susceptibility 
becomes ^ 

o \£f\ 

At this stage, we have different comments. First, it is easy to see that if C is constant then (jlUOp 
monotonically decreases as |ei?| increases. Second, let us write (jlOOp as 

which tells us that the first term is due to the disorder and the second is a manifestation of the 
confinement. Clearly, this can be interpreted as a quantum correction to the first term. According to 
pOip . this conclusion disappears if k = 0. Furthermore, in the limit C — )■ 0, the susceptibility becomes 
a (5-function ^ 

Xier) ~ -^SisF). (102) 
Obviously, this conclusion can not be reached in the clean limit Uj — t- as in |13j . 

5.3 Long-ranged disorder 

It is worthwhile to ask about the susceptibility for the long-ranged disorder. In a similar way to the 
short-ranged case, we use ()67p and ()72p to obtain 



X = --^ J ^defie) j Je' lm-^ [S+t(e', i3) + S,-,(e', i?)] 



(103) 

B=0 
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By analogy to ()75p . we change the variable to 



=e± T.i 



and define S^tt = iyti{X+ , , B) as 



Sit -^M EE 



n=0 fc=0 



St-ot-cMEE 



n=l A:=0 

The derivatives of Stot can be written in terms of Stot as 



X g{£n,k) 

X+X- - {en,kV 
X+X- -{en,k)^' 



dB 
dB^ 



dB 



52 



dB 



tot 2 tot "^^tot I ^ -^tot t^^tot ^-^tot 



dX^^dB dB dX^dX^ dB dB 



(104) 

(105) 
(106) 

(107) 
(108) 



where i, j, I = it and repeated indices indicate summation. The involved matrix elements are given by 

-1 

(109) 



dXi 



One can calculate the derivatives of at S = 0, in a similar way to the short-range case, and 
then obtain those for using (|107p and (jlOSp . to get 



r(S+ + S- 



dB' 



(l + a + 2/3)-^ 













( c 


X3 


\ 6 



+ 



1 



2/3 C"' 



where X = limB_j.o X+ = lims_s.o X and the two parameters are 

9{t) 



a 



f3 



2C 



tdt 



(•00 

2C I tdt 



X^-t^ 



lo {x^-t^y 

We have now derived all ingredients to write the expression for the susceptibility in this case: 



(110) 
(111) 

(112) 
(113) 



X 



defie) 



X(e) 1 

dX' lm-3 

X(-oo) ^ 



-l + SC 



I -a' (l-a')V. 



(114) 



B=0 



where the integration over e' has been replaced by 

dX' = [l + a' + 2^')'^ de'. (115) 
For the long-ranged disorder in the region |e| ^ (|112p and pi3p can be approximated as 



a 



/3 



C. 



(116) 
(117) 
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By substituting them in ()114p and performing the integration over X' , we see that the susceptibiUty 
is given by 



X 



Qttk 



def{e) Im 



1 

X2 



1 



3C 



1 + Clog 



IT' 



(118) 



B=0 



with X = e — Stot(£)- This can be rewritten as 



X 



/oo 1 2 foo 



X2 



1 + Clog 



(119) 



B=0 



The first term is the short-range contribution ()87p . The second term can be regarded as a contribution 
of order of 0{C), but this gives a minor effect since C is assumed to be small. When 0{C)'^ is 
neglected, the susceptibility becomes just 1 — 3C times as large as in the case of short-ranged disorder. 
Accordingly the integration of x over e weakly depends on C, while in C — t- we again get (jl02p . 
Finally, note that taking k = we end up with the results obtained by Koshino and Ando in |13j . 



6 Conclusion 

The present paper was devoted to give a complete solution to the confined Dirac fermion system in 
the presence of a perpendicular magnetic field. Using a similarity transformation we have formulated 
our problem in terms of the polar coordinate representation that allows us to handle easily the exact 
relationship between spinor components. One spinor component was obtained by solving a second 
order differential equation while the other component was obtained using the exact relationship ()2ip . 
It resulted in a complete solution space made of 8 subspaces, which suggests that it is necessary to 
include all components of this subspace in the computations of any physical quantity. A failure to do 
so will result in erroneous conclusions. 

Considering the disordered graphene, we formulated our solutions to capture the impurity po- 
tential. After providing the necessary tools, we used the Green function technique to determine the 
self-energy. For this, two cases are discussed, which concern short and long-ranged disorder. These 
allowed us to give a simplified form of the density of states that was used to determine different 
thermodynamical quantities. 

For further studies, we distinguished between two cases. As far as the short-ranged disorder is 
concerned, we further simplified the self-energy and found a quantum correction to the susceptibility 
for Dirac fermions in disordered graphene. On the other hand, we noted that «; = allowed us to 
recover the results of disordered graphene in the presence of a magnetic field [13]. Furthermore, in 
the limit C — > 0, the susceptibility becomes a (^-function (|lU2p . 

Subsequently, we treated the long-ranged case and stressed what makes a difference with respect 
to the former one. By comparing the results obtained in (|87p and ()118p . we noticed there is an extra 
term of the order of 0{C), but this gives a minor effect since C is assumed to be small. When 0{C)'^ 
is neglected, the susceptibility becomes just 1 — 3C times as large as in the short-ranged disorder. 
Accordingly the integration of x over e weakly depends on C, while in the limit C — )• we again get 
([102]). 
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The present work can be extended to other cases. One may numerically check the above results 
and compare them with those obtained before. On the other hand, an interesting question arises 
about what happens for a variable magnetic field, in particular an exponential variation of the field. 
This question is under investigation. 
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